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Abstract. 

The temperature T and magnetic field H dependencies of the electrical resistivity 
p of the recently discovered heavy fermion superconductor PrOs4Sbi2 have features 
that are associated with the splitting of the Pr 3+ Hund's rule multiplet by the 
crystalline electric field (CEF). These features are apparently due to magnetic 
exchange and aspherical Coulomb scattering from the thermally populated CEF- 
split Pr 3+ energy levels. The p(T) data in zero magnetic field can be described 
well by calculations based on CEF theory for various ratios of magnetic exchange 
and aspherical Coulomb scattering, and yield CEF parameters that are qualitatively 
consistent with those previously derived from magnetic susceptibility, specific heat, 
and inelastic neutron scattering measurements. Calculated p{H) isotherms for a T3 
ground state qualitatively account for the 'dome-shaped' feature in the measured p{H) 
isotherms. 



PACS numbers: 71.20.Eh, 71.27.+a, 74.70.Tx, 75.30.Mb 
1. Introduction 

The filled skutterudite compound PrOs4Sbi2 was recently discovered to be the first Pr- 
based heavy fermion superconductor, with a superconducting transition temperature 
T c = 1.85 K and an effective mass m* « 50 m e , where m e is the mass of the free 
electron [H Features observed in the magnetic susceptibility x(T), specific heat 
C(T), and inelastic neutron scattering (INS) |Hj all indicate that the Hund's rule J = 4 
multiplet of the Pr 3+ ion is split by the crystalline electric field (CEF). An analysis of 
these data within the context of a cubic CEF using the theory of Lea, Leask and Wolf 
(LLW) yielded a T3 nonmagnetic doublet ground state, a T 5 triplet first excited state 
at ~ 10 K, and higher T 4 triplet and Ti singlet excited states at ~ 130 K and ~ 313 K, 
respectively. A 'roll-off' near ~ 8 K with electrical resistivity p(T) in zero magnetic field 
was attributed to a decrease in spin or charge dependent scattering from the low-lying 
T 5 excited state whose population decreases with decreasing temperature, while features 
in p(H,T) at low temperatures < 1.5 K and high fields > 4.5 T were identified with a 
high field ordered phase (HFOP) of magnetic or quadrupolar character [3 El- 



Crystalline electric field effects in the electrical resistivity of PrOs^Sb^ 



2 



In a lanthanide compound, the conduction electrons interact with the localized 
4f electrons through magnetic exchange and the direct Coulomb interaction. For 
lanthanide ions whose 4f electrons are hybridized with conduction electrons, the 
exchange interaction can be negative, which favors screening of the magnetic moments 
of the 4f electrons by the conduction electron spins. This leads to contributions to 
the electrical resistivity that increase as log T with decreasing temperature that are 
associated with the Kondo effect. Since no Kondo-like contribution has been found 
in p(T), we disregard 4f conduction electron hybridization and treat the 4f electrons 
in the ionic limit, as a first approximation. When the degeneracy of the Hund's rule 
multiplet of a lanthanide ion is lifted in a crystalline electric field, the change of the 
resulting lanthanide energy level populations with temperature introduces temperature- 
dependent features in the electrical resistivity. These features are more subtle than 
those seen in properties such as magnetic susceptibility and specific heat, but can still 
be described based on the CEF Hamiltonian of Lea, Leask, and Wolf (LLW) [5 . 

This paper presents calculations of the electrical resistivity due to magnetic 
exchange and aspherical Coulomb scattering associated with the CEF splitting of the 
Pr 3+ Hund's rule multiplet in PrOs4Sbi2- Electrical resistivity p(T) data in zero 
magnetic field are directly fit with the calculated equations, providing parameters that 
are used to calculate p(H) isotherms. It is shown that a T 3 ground state can qualitatively 
account for the 'dome-shaped' features that have been observed in the measured p{H) 
isotherms. 



2. Calculation of the electrical resistivity in the crystalline electric field 

The effects on electrical resistivity due to magnetic exchange [HI IZ] and aspherical 
Coulomb |H| scattering have been separately considered many times in the past; often 
one effect was neglected over the other. The present calculations follow the procedure of 
Fisk and Johnston [9], where it was demonstrated that both contributions are important 
for describing the temperature dependence of the electrical resistivity for the compound 
PrBg and, in addition, the relevant equations were presented in a form appropriate to the 
Stevens operator equivalent formulation of the LLW crystal field Hamiltonian [TUl lllj . 
In this notation [7j, the total contribution to the electrical resistivity from CEF effects, 
Pcef, is: 

Pcef = Po[rTr(PQ M ) + (1 - r)Tr(PQ A )], (1) 

where r is a coefficient representing the ratio of the magnetic exchange term to the 
aspherical Coulomb scattering term. The temperature dependent matrix, Py, is common 
to both terms: 

p. e~™ m - Ej) 

where Ei are the eigenvalues of the CEF eigenstates and (3 = 1/ksT. The 
matrix represents magnetic exchange scattering, and the Qf^ matrix is associated with 
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aspherical Coulomb scattering due to the quadrupolar charge distribution of the Pr 3+ 
ion. The matrices are: 

^ = I^V,li)r + ^l(^ + li)l 2 + ^l(^-li)l 2 , (3) 

Qt 3 = E KW>l a . (4) 

m=-2 

where the |i)'s are the CEF eigenstates, and the |/™'s are the operator equivalents of the 
spherical harmonics for L = 2 (i.e., quadrupolar terms), and are given elsewhere |12j . 
The Qij matrices are also normalized to each other [Oj, such that 

£<$ = J2Qtj = (2J + 1)J(J+ 1) = 180 (for J = 4). (5) 

Using Matthieson's rule, the electrical resistivity is separated into impurity, phonon, 
and CEF contributions: 

P = Pimp + Ap La + pcef (6) 

where p imp is the impurity scattering term, p La is the lattice term of the previously 
measured isostructural compound without f-electrons, LaOs4Sbi2 [H|, and pcef is 
given by Eq.fpQ). The constant A was used to scale the lattice contributions of p(T) 
of LaOs4Sbi2 and PrOs4Sbi2 to one another, assuming they have the same temperature 
dependence, and to account for uncertainties in the geometrical factor due to the 
irregular shape of the crystals and/or microcracks in the crystals. For Pcef, the energies 
Ei can be expressed in terms of the parameters x and W in accordance with the LLW 
formalism [Sj. In the absence of a magnetic field, the eigenstates \i) are the same for 
all values of x and W for J = 4 in a cubic CEF. Additional terms due to tetrahedral 
symmetry in the crystal field Hamiltonian ^3] were not included in this analysis in 
order to keep the calculations consistent with those previously made on these samples. 
Eq. (JSJl was fit to the original data from Ref. [2] between 1.9 K and 20 K, assuming a 
T 3 doublet ground state, with results shown in Fig. [JJ The best fit to the data yields an 
impurity scattering term p imp = 1.67 pfl cm, a scaling factor for the LaOs 4 Sbi 2 lattice 
A = 0.21, a CEF scaling factor po — 0.385 pQ cm, a CEF scattering mechanism ratio 
r = 0.46, and LLW values of x = —0.7225 and W = —2.97, which results in the first 
excited state T 5 triplet lying ~ 5 K above the T 3 doublet ground state. 

It is notable that the fit of Eq. |U1 to the zero field p(T) data between T c and 
20 K is excellent (upper inset of Fig. EJ, and the agreement between the fit and the 
data at high temperatures > 50 K is qualitatively good (lower inset of Fig. [TJ. Since 
the magnetic exchange and aspherical Coulomb scattering contributions to the CEF 
are nearly temperature independent above ~ 75 K, the discrepancy is most likely due 
to the LaOs 4 Sbi2 lattice term not being an accurate representation of the PrOs 4 Sbi 2 
lattice term, which is reasonable considering the anomalous resistivities of other La- 
based compounds compared to their Y and Lu-based counterparts without f-electrons 
[To] . Unfortunately, single crystals of these YOs 4 Sbi 2 or LuOs 4 Sbi 2 compounds, which 
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Figure 1. Electrical resistivity p vs temperature T between 1 K and 50 K for a single 
crystal of PrOs4Sb 12 . The solid line is a fit of the function described in the text (Eq.|SJ, 
consisting of contributions from the scattering of electrons by impurities, phonons 
(represented by the scaled phonon contribution to the resistivity of LaOs4Sb 12 ), and 
magnetic exchange and aspherical Coulomb interactions. Upper inset: p vs T between 
1 K and 20 K for the same single crystal of PrOs4Sbi2. The dashed line is the CEF 
contribution to the fit represented by the solid line. Lower inset: p(T) of PrOs4Sbi2 
between 1 K and 300 K. The CEF contribution (dashed line) is nearly negligible at 
higher temperatures. 

may provide more accurate MOS4SD12 filled skutterudite lattice resistivities, are not 
available. The data were also fitted with the CEF scattering mechanism ratio r fixed 
at intervals of 0.25, and the results are presented in Table ^ All values of r resulted 
in satisfactory fits of the p(T) data at low temperatures, with the value of the splitting 
between T 3 and T 5 , A 35 , increasing from 3.3 K for r = 1 (magnetic exchange scattering 
only) to 6 K for r = (aspherical Coulomb scattering only). The values for r = 
are closer to the values determined from x{T), C(T), and INS measurements. It is 
possible that a better determination of the lattice term could result in a more definitive 
result for the value of r. While it is unlikely that magnetic exchange scattering is 
completely unimportant, it does appear that aspherical Coulomb scattering contributes 
considerably to the electrical resistivity of PrOs 4 Sbi 2 . 
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Table 1. Parameters resulting from a fit of Eq. El to the data in Fig. ^ The value 
of the CEF scattering ratio r was fixed at various values to determine how the other 
parameters would change. The larger the aspherical Coulomb scattering contribution 
(smaller r), the closer the value of A 35 to values previously calculated. 



r 


Pimp 


A 


Po 


X 


W 


A35 




(fifl cm) 




(fitt cm) 




(K) 


(K) 





1.88 


0.204 


0.373 


-0.716 


-2.52 


5.99 


0.25 


1.80 


0.207 


0.378 


-0.720 


-2.78 


5.33 


0.5 


1.64 


0.208 


0.387 


-0.723 


-3.00 


4.73 


0.75 


1.33 


0.206 


0.406 


-0.726 


-3.18 


3.93 


1 


0.86 


0.200 


0.435 


-0.728 


-3.30 


3.33 



3. The high field ordered phase in PrOs 4 Sb 12 

Measurements of the electrical resistivity, specific heat, and thermal expansion of 
PrOs4Sbi2 at low temperatures < 1.5 K and high magnetic fields > 4.5 T reveal features 
that appear to be associated with an ordered phase, one that may be quadrupolar in 
nature [3J EJ H3 HE] • In electrical resistivity, these features are manifested as 'roll- 
offs' in p(T) and large 'domes' in p(H) (31 HH1 UH] • The 'domes' in p(H) are especially 
intriguing, as it is expected that scattering should be lower in an ordered state. It has 
been suggested that the crossing of the CEF energy levels in magnetic fields are related 
to the high field ordered phase [El HHj and can contribute to the increase of scattering 
in this region. Recent neutron diffraction experiments have also suggested the existence 
of an antiferro-quadrupolar ordered state that is associated with magnetic field-induced 
energy level crossings for a Ti ground state [2D]- To check this hypothesis for electrical 
resistivity, p was calculated as a function of H in order to qualitatively observe the 
shape of the CEF contribution to the electrical resistivity. 

Shown in Fig. a) are the energy level splittings as a function of magnetic field 
(H || (001)) for the T 3 doublet ground state and T 5 triplet first excited state, using the 
LLW values of x and W derived from the fit to the electrical resistivity of PrOs4Sbi2 
with the CEF scattering mechanism ratio r chosen to be 0.5. This value of r was used 
because it is consistent with the value determined by Fisk and Johnston for PrBg [0] 
as well as being a compromise between aspherical Coulomb scattering, which results 
in a A 35 that is closer to values estimated from other analyses (x(T),C(T), INS), and 
magnetic exchange scattering, which is more universal in lanthanide compounds. The 
corresponding plot of p vs H for various temperatures is shown in Fig. E^b). The 
parameter po was chosen so the calculated p{H) curves were on the same scale as 
those measured experimentally. A dome shape, which sharpens considerably at lower 
temperatures, can clearly be seen in the graph. A plot of p vs H for PrOs4Sbi2 at 1.4 K 
[To] is presented for comparison to the theoretical curves. The maximum of the dome 
corresponds to the the magnetic field at which the ground state changes from one of 
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Figure 2. (a) Low temperature energy level splittings £ as a function of magnetic 
field H for the LLW parameters x = —0.723 and W = —3.00. (b) Calculated electrical 
resistivity p vs magnetic field H for the energy splittings in (a). It was assumed that 
magnetic exchange and aspherical Coulomb scattering contributed equally to the CEF 
resistivity (r = 0.5 in Eq. QJ. The calculated lines were also scaled by po = 0.5 in 
order to bring the curves to the same scale as those measured experimentally. For 
comparison, the dashed line is the measured p vs H for PrOs4Sbi2 at 1.4 K Jfaj- ( c ) 
E vs H for the LLW parameters x — 0.508 and W — 2.46. (d) p vs H for the energy 
splittings in (c), and similar constraints as (b). 



the T 3 states to one of the T 5 states. The dome structure can be understood by taking 
the limit of the equation for the terms, Eq. as the lowest excited state Ei — > 0, 
assuming Ej = is the ground state, and the other energy levels are large enough to be 
ignored in the sum in the denominator. Eq. |2] then becomes 

The shape of this equation as a function of energy is a dome centered on E$ = that 
sharpens considerably as the temperature is lowered. Since Ei is nearly proportional to 
H (Fig. E(a)), domes also appear in p(H). These domes are present in both magnetic 
exchange and aspherical Coulomb scattering, as there are T 3 — > T 5 transitions in both 
mechanisms, which unfortunately contributes to uncertainty in determining the relative 
importance between the two types of scattering. 

Similar to CEF fits to magnetic susceptibility x{T), it is also possible to fit p(T) 
for PrOs4Sb 12 at low-temperatures T < 20 K with an energy level scheme that has a 
T 1 singlet ground state and a T 5 triplet first excited state. The fit resulted in similar 
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impurity and lattice fitting parameters as for a T 3 ground state, with LLW values of 
x = 0.508 and W = 2.46. The parameters from the zero-temperature fit were used to 
calculate the energy level splittings as a function of magnetic field, as shown in Fig.EJc), 
with the corresponding plot of p vs H shown in Fig. |2Jd). It is immediately apparent 
that the dome shapes in p(H) due to CEF scattering are less pronounced and are not in 
agreement with the data when the ground state is Ti, even though the fits to p(T, H = 0) 
are nearly identical. A qualitative reason for this is that Ti — > T 5 transitions only occur 
in aspherical Coulomb scattering, and the energy levels cross at a much lower field than 
for the T 3 ground state situation. 

4. Summary 

Features in the electrical resistivity of PrOs4Sbi2 can be explained in terms of effects 
arising from the temperature- and magnetic field-dependent populations of crystalline 
electric field energy levels. Fits to the p(T) data were accomplished by a combination 
of magnetic exchange and aspherical Coulomb scattering in the context of a r 3 ground 
state. The parameters resulting from this fit were used to calculate isotherms of p{H) 
that agree qualitatively well with isotherms previously measured. While it was possible 
to fit the p{T) data well with a Pr 3+ energy level scheme with a Ti ground state, when 
p(H) isotherms were calculated using CEF parameters for this case, they were not in 
good agreement with the measured p{H) isotherms. It is evident that the high field 
ordered phase that is also observed in high-field measurements of C(T) and a(T) is 
closely related to the ground state crossover in the energy level-magnetic field phase 
diagram. 
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